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Abstract 

The structures of order parameters which determine the bounds of the phase 
states in the framework of the CP^ Ginzburg-Landau model were considered. Using 
the formulation of this model in terms of the gauged order parameters (the unit 
vector n, density p and momentum of particles c) we found that some universal 
properties of phases and field configurations are determined by the Hopf invariant, 
Q and its generalizations. At a sufficiently high level of doping it was found that 
beyond the superconducting phase the charge distributions in the form of loops 
may be more preferable than those in the form of stripes. It was shown that in the 
phase with its mutual linking number L < Q the transition to an inhomogeneous 
superconducting state with non-zero total momentum of pairs takes place. The 
universal mechanism of the topological coherence breaking of the superconducting 
state due to a decrease of the charge density was discussed. 

PACS 74.25.-q, 74.80.-g, 71.10.Hf, 71.10.-w 

1. Introduction 

Among the challenging problems of cooperative phenomena in planar systems near Mott 
transition, there are those which at the first glance may not be associated with the 
appearance of high-temperature superconducting states in doped antiferromagnetic insu- 
lators. For example, we are interested in the origins of qualitatively similar cooperative 
behavior in various compounds and very rich content of their phase diagram, as well 
as in origin of the emergence of inhomogeneous states, the typical for such systems 
[2l IHl Bl ID Low-dimensional structures in the distribution of spin j2lE| and charge 
00 El degrees of freedom exist in the state which precedes the high-temperature super- 
conducting phase. Keeping in mind this property [3 we should choose such a model for 
describing the mentioned phases which will contain them as limiting cases. The mean 
field Ginzburg-Landau theory with the appropriate choice of the order parameters may 
be used for understanding of such a kind general problems. The key question in this 
universal approach is the method with the aid of which the order parameters encode 
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simultaneously the content and the distribution of charge and spin degrees of freedom 
of excitations in various phase states. 

The recent progress in solving analogous problems in the non-Abelian field theory 
[H] and its development in condensed matter physics has shown that CP^ Ginzburg- 
Landau model is preferable. The two-component order parameter of this model is used 
for solving the problems of two-gap superconductivity P] . In the theory of electro- weak 
interaction it has the sense of the Higgs doublet of the standard model. In this paper 
we will assume that the order parameter is a spinor realizing two-dimensional representa- 
tion of the braid group which arises due to classifying the quantum states at permutations 
of particles in (2+l)-dimensional systems. Considering factorization with respect to the 
centre of this non-Abelian group we obtain the gauged CP^ Ginzburg-Landau model. 
We pay attention that the order parameter of this model is two-dimensional 0^]. Only 
in this case we can introduce the unit vector field which describes the distribution of one- 
half spin degrees of freedom in the long-wavelength limit as well as use the Hopf invariant 
for classifying n-field configurations and consider correctly the phases with different dis- 
tribution of charge degrees of freedom. Due to the above mentioned reasons we use the 
generalized n-field model which after the exact mapping of CP^ Ginzburg-Landau model 
^ includes Faddeev term pJI- Because of the non-Abelian gauge theory origin of this 
significant part of the model we hope that the obtained answers are universal and will 
give a deeper insight into the discussed problems. 

The Hopf invariant describes the degree of linking or knotting of the filamental man- 
ifolds where the field of the unit vector n is defined. The study of the behavior of the 
vortex filament tangle is a separate problem and attracts attention due to several rea- 
sons. First, at small distances the topological order associated with the linking exists 
against the background of the disorder caused by arbitrary motion of separate parts of 
the system of entangled vortex filaments. Thus, unlike point particles, the properties of 
the tangle are determined by the behavior of its fragments in the ultraviolet and infrared 
limits. Since the coordinates of the vortex core are canonically conjugated, the mentioned 
circumstances are presented by noncommutativity of these variables depending on the 
linking degree. Second, such soft medium as the tangle of linked filaments is the hot 
problem in condensed matter physics and beyond, in particular, in connection with the 
DNA problem. After the coupling constants are transformed into the searched functions 
the appearance of the double helix as the solution of the equations of motion in the soft 
variant [E] of the n-field model is its general property. 

In the present paper we consider some properties of field configurations in CP^ 
Ginzburg-Landau model defined in the following paragraph. The main goal of the paper 
is to find the bounds of free energy in the superconducting state and in the inhomoge- 
neous phase with broken antiferromagnetic order, as well as to describe the properties 
of the charge density distributions corresponding to this state. Considering the non- 
superconducting phase in the soft version of the model we analyse the value of the 
contribution to the free energy of the charge density distributions in the form of loops 
and stripes. Along with the results from brief publications |14[ 1 151 IT?!] in the third and 
fourth sections, we discuss the properties of the inhomogeneous superconducting state 
with nonzero total momentum of pairs and compare it with the LOFF states ^3 E| 
and with the results from the recently proposed BCS-like model with two types of 
particles. We also pay attention to the dependence of the bounds of phase states on the 
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generalized (2+l)D Hopf invariant in the case of S*^ x 5*^ — > 5*^ and x x ^ 
mapping classes and on the external magnetic field. In conclusion we discuss some open 
problems. In appendix there is the proof of the inequality which determines the relation 
between the contributions to the free energy of n- and c-field configuations. 



2. CP^ Ginzburg-Landau model 



We will use Ginzburg-Landau model: 
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with a two-component order parameter, 
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(1) 



(2) 



which satisfies the condition of + |x2p = 1- This constrain of two components 
Xa takes place in the complex projective space CP^, for which the considered model 
is defined. The model (1), (2) with different masses was used before E] [1,6] in 
the contex of two-gap superconductivity, as well as in the standard model of the non- 
AbeHan field theory [HlEl- In the present paper we consider the states in planar systems. 
Thus we suppose that ^ has the sense of the order parameter reaHzing two-dimensional 
non-Abelian representations of the braid group used for classifying quantum states at 
permutations of particles in systems with two spatial dimensions. Realizing Abelian 
projection [201, the vector Ak compensates the local choice of the phase of the function 
The last terms in (1) describe Ginzburg-Landau potential y(\E'i,^'2) and the self- 
energy of the gauge field. It has been shown recently [T] that there is an exact mapping 
of the model (1), (2) into the following version of the n-field model: 
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The free energy in Eq. (3) is defined by the scalar - the density of particles , the field 
of the unit vector = xcr^x, where x = {XiiXi)^ is Pauli matrix, and the field of 
the momentum c = J/p^ = 2(j — 4A). The total current J contains paramagnetic part 
j = i[xi'^Xi ~ c-c. + (1 ^ 2)] and diamagnetic term — 4A. Equation (3) was obtained 
with the use of the following notations: Fik = diCk — dkCi, Hik = n • [din x i9in], and 
dimensionless units as follows: of the length L = (fi -l-^2)/2 with the coherence length 

— h/^J2m ba, the momentum h/L as the unit of the momentum c, the density units 
c^/(5127re^L^) of particles per mass unit (at parametrization of '^a in the form (2)) and 
the energy units j/L with 7 = {ch/ef /5127r. 

In formulation (3) the Ginzburg-Landau functional depends on gauged order param- 
eters p^, c, and n. They characterize spatial distributions of charge and spin degrees 
of freedom with current or without it. The functions Xa determine the orientation of 
the unit vector n which describes (in the long-wavelength limit) the properties of the 
magnetic order. Besides, functions Xa define the value of the paramagnetic part of the 
current. Comparing different forms of presentation of CP^ Ginzburg-Landau model we 
note that vortex field configurations ^'q in the model (1), (2) are equivalent to textures 
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of the field n in terms of the model (3) . We also pay attention to the fact that the ansatz 
(2) has the sense of factorization of longitudinal p and transversal Xa degrees of freedom. 
In the superconducting state the composition of spin j and charge degrees of freedom is 
important, since the current contains diamagnetic U{1) gauge component — 4A. 

In the soft variant of the extended model of n-field (3) the multipHers of the first 
term describe the distributions of spin stiffness and the square of the inverse length 
of the density screaning. It is seen from this example that the competition of the order 
parameters p, n and c may be the origin of the existence of the phase states with different 
ordering of charge and spin degrees of freedom. We enumerate the limiting cases of the 
model (3) in inhomogeneous (n ^ const) situations: 

1. A state with a broken antiferromagnetic order: c = 0, p = const. 

2. A state with a quasi-one-dimensional density distributions: c ^ 0, p ^ const. 

3. An inhomogeneous superconducting state: c ^ 0, p ^ const. 

4. c ^ 0, p 7^ const. 

In the case of n = const and c ^ 0, p ^ const functional (1) is equivalent to the 
one-component Ginzburg-Landau model. 



3. The bounds of the free energy 

1. A phase state with a broken antiferromagnetic order. - Let us consider the first case 
in the mentioned list. In this limit the free energy equals 

j [51 {dunf + 52 (n • [d,n x S^n])'] . (4) 

We supposed that in the considered phase the constant value p — po may be found from 
the minimum of the potential V and introduced the notation gi for the coupling constants. 
The properties of the model (4) were studied in detail in Refs. |2T1 l22l l23l |24. ,25 . .26. .27) . 
The analysis of the dimensionality shows that the first term in Eq. (4) is proportional to 
the characteristic size Rq of the n-field configurations, and the second term is inversely 
proportional to this scale. Therefore, the energy (4) has a minimum which is achieved at 
Rq = Vys/ffi- This explains why the second term in Faddeev-Niemi model (4) allows 
to avoid Derrik's restriction of the existence of three-dimensional static configurations 
with the finite size. In the infrared limit this term characterizes the mean degree of 
noncollinearity < 0|Si • [S2 x S3] |0 > in the orientation of three spins which locate the 
sites of the quadratic placket. 

It was shown in Refs. |25[ I2(il l?7] that the lower energy bound in the model (4) 

F ^ 32^2 |g|3/4 (g) 

is determined by the Hopf invariant 

Q = J d^xEikiaidkai . (6) 

In this equation the vector denotes the gauge potentrial which parametrizes the mean 
degree of the noncoUinearity of the orientation of neighboring spins in the following way: 
Hik =n - [diU X dkn] = d^Uk - dkUi. 
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The dependence Fmin IQP^'* (5) with the boundary conditions n (0,0, 1) in the 
space infinity was verified in Refs. |22l 1231 124j in simulation of the configurations of 
n-field. Such a boundary condition means that the space of the n-field definition 
effectively compactizes into a three-dimensional sphere Thus the unit vector n 
accomplishes the mapping of S'^-sphere into the space of the two-dimensional sphere §^ . 
Let the vector n be directed to some general point of a two-dimensional sphere. We will 
be interested in the answer to the following question: what is a pullback of this point 
in the space Or in other words, what set of points from the definition field of the 
vector n{x,y,z) contributes to the point of the target two-dimensional sphere? Since 
the space §^ is compact and its dimensionality is greater by unity than that of the 
sphere §^ , the pullback of points on the sphere §^ are closed and in general linked lines 
on the sphere S'^ . The Hopf invariant Q (6) describes the degree of linking or knotting 
of these lines. It belongs to a set of integers Z to which the considered homotopic group 
7r3(S^) = Z equals. In particular, for two once-linked circles Q — I equals 1, for one of 
the simplest knots (a trefoil) Q equals 6 and etc. As a result, the n-field configurations 
are divided into classes corresponding to the values of Hopf invariant. We should 
emphasize once again that linked or knotted configurations may be numbered by the 
Hopf index only in the case of the CP^ Ginzburg-Landau model with its two-component 
order parameter, because at npn M > 1 the homotopic group ^^{CP^^) = is trivial . 

2. Quasi- one- dimensional density distributions. - Let us consider the states outside 
the superconducting phase from the second line of the list of limiting cases, to which the 
CP^ Ginzburg-Landau model leads. In this soft version of the model the functional (3) 
has the form 
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(7) 



In Eq.(7) the positive constant b correponds to the phase with the broken antiferromag- 
netic order. 

The state with the broken antiferromagnetic order considered above has a lower 
energy than the "soft" state we are interested in now. The latter may be metastable 
[28| . In this section, we will consider just such states and compare their contribution 
to the Ginzburg-Landau energy without studying the problems of their relaxation, the 
critical sizes of nuclei of different phases and etc., which are of separate interest. 

Under the condition that the electron spin and charge are transferred from one of four 
sites of some plaquettes to the dopant reservoir, the terms with Hik in Eq.(l) characterize 
(in the infrared limit) the mean degree of non-collinearity < 0|Si • [S2 x S3] |0 > in the 
orientation of three spins, which remain in the sites of a quadratic lattice plaquettes. 
Therewith the deficit of the charge density relates to the density p^, describing in 
(1) the distribution of the exchange integral, by the relation + p\= const. From the 
long- wavelength point of view the distribution of the spin density p^ in the limited region 
with the exponential law of decrease at its boundary (for example, at the distribution 
in the circle with radius tq with the exponential decrease at the length R <^ tq) will be 
accompanied by the quasi-one-dimensional distribution of the charge density p\ along 
the boundary of this region, i.e. along the ring with thickness R and radius ro. It is seen 
from here that studying spatial configurations of the density field p^ in planar systems 
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makes it possible to find the form of one-dimensional distributions of the electric charge 
density with the aid of the above-mentioned holographic projection. 

It has been known for a long time that in such a phase state the distributions of 
the charge density {dup)^ have the form of stripes^. Due to the gradient term {dkpf 
in (7), quasi-one-dimensional field p configurations are really preferable. It seems to 
be almost obvious also that density distribution in the form of rings give the smallest 
contribution to the energy. Let us find the contribution to the free energy (7) from quasi- 
one-dimensional density distributions p^ in the form of rings and stripes and compare 
the computation results with the experimental data. We will choose the following trial 
functions for the field p configurations in the form of a ring and a stripe: 

p = poexp[-{r-rof/2R^] (8) 

and 



Here po = yb/d, tq is the ring radius, R is its width, 2Ly = 27rro is the stripe length, 
Lx = R is its width. Since configurations (8) and (9) do not depend on the third 
coordinate, we will assume that along it the size is limited by the length and also that 
i? < To. The calculation of energy (7) with the aid of (8) and (9) yields the following 
results for the contribution to the free energy from the ring Fr (at R <^ ro) and the 
stripe Fxy-. 



^pIl..4[i + ^] , (10) 



Here fo = V^ro, 1/^ = 2 [no - (1 - 1/^8)6] , no is a certain characteristic value of the 
"multiplier" (Ofen)^ in (2), which is of the order CiR~^, whereas b = C2R~^ST, where 
c, ^ 1 and ST = [Tc — T)/Tc. In these equations we omitted the term Hf^. from Eq.(7) 
since we consider that it is approximately the same for the both types of distributions. 
The exact equation for Fr (in units ttpqLz) contains the term 
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5Fr = 2[h{xo) - xol2{xo)\+'^ h{xos/2) - xos/2Io{xos/2) -h^ [/i(a;o) - a;o7o(a;o)] 
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But it is exponentially small already at R/tq ~ 1/4 and i? ~ A with 1/A^ = 2(no — b): 
SFr ~ 10-^. 

For the optimal width i? = ^ (at i? <C tq) the difference in free energies AF = F^y—Fr 
in units np^Lz has the form: AF = 1 -|- ci — (3/4)c2(5T. One can see from this 
equation that at 4(1 + Ci)/3c2 < 1 in the temperature range [1 — 4(1 + Ci)/3c2] Tc < 
T < Tc bordering on the critical temperature Tc of the transition to the state with 
the spin pseudo-gap, rings are preferable (see Fig. 1). In the temperature range 
T < Tc [1 — 4(1 -I- ci)/3c2] stripes are the main configurations. As it is known one 
may approach T,. keeping the temperature constant by increasing the level of doping. 

^We suppose that the stripe characteristic sizes are essentially greater than the lattice scale. In this 
case, the use of phenomenological approach of the Ginzburg-Landau mean field theory is justified. 
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Figure 1: Schematic representation of closed (b) and open (c) quasi-one-dimensional 
structures of the charge density (see. j22]) around antiferromagnetic dielectric nano- 
clusters (a). 



The recent paper [5^1 gave the evidence of the existence of ring-shaped charge struc- 
tures obtained just under such experimental conditions. In a certain sense, the tunnel 
microscope in this experiment (221 cohects data from a two-dimensional shoe of knots 

m 

Let us make several remarks concerning the spin density distribution in the disk, 
surrounded by a ring charge distribution. The above-mentioned statement concerned the 
case when spin disorder arose only in the region directly adjacent to the disk edge and, 
as a result, we had an antiferromagnetic phase inside. If the antiferromagnetic order is 
broken everywhere in the disk, it is necessary to consider the corresponding distribution 
of the density in the form of a disk in order to compare its contribution with Fr. 
When we consider the contribution of the density distributions to the free energy in 
the form of a disk, we get a double gain (in comparison with rings) due to the existence 
of one edge and have a loss due to the area. The calculation shows that at small R/vq 
the distributions in the form of rings appear to be more preferable. 

Let us consider the dependence of the critical temperature Tc on the level of doping. 
To do this, we present the relation of Fr to Fxy in the following form: 

(12) 



where 

2 



Fxy l + BR/fa' 

3 710 - 0.68 5 



^^1 no -3/46 



4 no -0.65 6' 



2 no - (1 - l/V^) b\ 

One can see that configurations in the form of stripes are more preferable in the 
range 0.65 < no/6 < 0.68, where F^y < Fr. Normalizing the density Pq to the 



particle number N by the condition N = 2Trro^Lzmpl, we get the relation between 
the parameters uq and b, which we write in the form of a; = b/^/no — 0.655. Here 
X — Nd/{^/2rmrLzro). Thus, for the bounds of the considered above range, where 
no ^ b ^ 6T, we have T{x) = Tc(l — Ax^) with a certain constant A. Therefore, inside 
the region, belonging to the phase state with the broken antiferromagnetic order, there 
is a narrower region, locating between the parabolas T{x), where the charge structures 
have the form of stripes. 



3. The inhomogeneous superconducting state. - We consider a superconducting state 
with finite value of the total current J which exists against the background of a certain 
n-field distribution, assuming that p ~ po = const. In this case the free energy equals 

F = Fn + Fc.- F„,t = (13) 



'3„ 



The negative sign of the interaction energy Fint of c- and n-fields appears because of 
diamagnetism of the considered state. As a result the coupling constant 172 = 1 of the 
term Hff. decreases due to renormalizing in such a way that energy of the superconducting 
state c 7^ is less than the minimum value in the inequality (5) . To find the exact lower 
bound of the free energy in the superconducting state c ^ 0, we will use the following 
auxiliary inequality 

(32^2)4/3 ^ (14) 

where the invariant 

i = J d^xeikiCidktti (15) 

determines the degree of the mutual linking [HOI 1^ of the current lines and the lines 
of the magnetic field H = [V x a] . Like Q it is the integral of motion |31l I32j in the 
considered barotropic state. The proof of inequality (14) [14 is given in the appendix of 
the paper. 

Linking indices, which characterize the correlations of spin and charge degrees of 
freedom, form the following matrix: 

= TTT^ / d'xe^kKdkaf = f ^, ) . (16) 



167r2 J ~ \ L Q' 

In this symmetric matrix (L = L') with aj = Ui and of = Ci the integral might be deter- 
mined also by the asymptotic Hnking number (201 ■ Let us pay attention to a circumstance 
which will be important below. Being normaHzed to the charge density, unHke the unit 
vector n, the vector of the momentum c = J/p^ belongs to the noncompact manifold. 
Because of this the Hopf numbers, defined with the aid of it in (16) are not integers in 
general: {L, Q') ^ Z. In the superconducting state, where Abelian U{1) gauge symmetry 
is broken and the charge is not conserved, the numbers L and Q' play the role of contin- 
uous interpolation parameters, which unite the considered compressed and incompressed 
{Kap € 1) phases. From this point of view, the superconducting states with Kap e Z 
and Kaf3 ^ Z belong to one and the same class of universality (HI . 
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To find the lower bound of the functional (13), alongside with Eq.(14) we will use 
Schwarz-Cauchy-Bunyakovski inequality: 

F,nt < 2||F,fe||2 • \\H,kh ^ '^Fl'^Fl'^ ■ (17) 

1 /2 

Here ||fifc||2 = [/ d?xFf^\ . We note that the equality in the r.h.s. of Eq.(17) is 
achieved in the ultraviolet limit, when the size of linked vortex configurations is small 
enough. Substituting the boundary value Fint in (13) we get 

F ^ Frmn = {Fl'^ - Fl'^f . (18) 

The Hopf configuration with (5 = 1, for which the lower limit takes place in Eq.(5), 
presents two linked rings with radius R and {Fn)min = 2'k'^R^{%/ + %/ R'^)\r^i — 327r^. 
We will assume that also in our case c 7^ there are configurations for which the equality 
in Eq.(14) is valid. Let us emphasize an important circumstance, which we will discuss 
more thoroughly in the next section. For small values of p and, therefore, for great values 
of the field c (since all terms in (13) are of the same order) we encounter the instability 
of linked configurations with respect to small perturbations. This leads to the restriction 
of values Fc from above. Keeping in mind this remark and using in Eq. (18) for Fc of 
the lower bound F^^^ = (327r2)4/3 i?-5/6 ^^ ^^^^ Eq.(14) and Fn = 327r^ \Q\^^'^, we get 
for the states with Q 7^ 0, that 

F ^ 327:' \Q\'/* {1 ~ \L\/\Q\f . (19) 

One can see from Eq.(19) that for all numbers L < Q the energy of the ground 
state is less than in the model (4), for which inequality (5) is valid. The origin of the 
energy decrease may be understood by comparing the values of different terms in Eq.(13). 
Even under the conditions of a considerable paramagnetic contribution j to the current, 
the diamagnetic interaction in the superconducting state for all classes of states with 
L < Q reduces in (13) its own energy Fc of the current and a part of the energy F„ 
associated with dynamics of n-field. In the considered state the total momentum of 
superconducting pairs c does not equal zero. In this respect the inhomogeneous state 
with current is analogous to the state proposed in Refs. [TTllTs]. 

4. The properties of phase states 

The phase state with the broken antiferromagnetic order at {dkp)' 7^ is a background at 
which the transition to the inhomogeneous superconducting phase with Fik ^ occurs. 
It is convenient to discuss the characteristics of this transition at the density change 
beginning from the superconducting state. In this phase the constant value of the charge 
density, due to the breaking of the gauge invariance C/(l), plays the role of the tuning 
parameter of the system. 

Let the parameter po change in some range. Since all terms in Eq.(3) are of the 
same order, as po increases, the momentum c and, consequently, the index of the mutual 
linking L decrease. In this case at sufficiently small L the smallest superconducting gap 
decreases with increase of Q against the background of a great value of 327r^ jQ^/'' of the 
spin pseudo-gap. 
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When po reduces, the following effect takes place. Being proportional to 5^ ^ Pa i 
the radius TZ of compactification M? S'^ grows till it exceeds some critical value TZcr ■ 
At TZ > Her the Hopf mapping is not stable relative to small perturbations of 
linked vortex field configurations. As a result, the symmetry U{2) which is associated 
with identical Hopf mapping appears to be spontaneously broken. This means that 
the topological configurations of field n and c, instead of being spread out over the 
whole space localize around a particular point (the base point of the stereographic 
projection R^^ ^ S'^) and collapse to locaHzed structures j^. We can see that there is 
an optimal value of po and, consequently, the values of the characteristic momentum c 
and the relation |L|/|(3|, for which there arises the greatest gain at the transition to the 
superconducting state. 

Till now the vector A has characterized the degrees of freedom, associated with the 
internal charge gauge symmetry U(\). If the external electromagnetic field is applied, 
the vector A equals the sum of internal and external gauge potentials. In the exter- 
nal magnetic field, due to diamagnetism of the superconducting state, the momentum 
c decreases. Like in the case of po increasing, this leads to suppressing the supercon- 
ducting gap. Playing the role of the smooth tuning parameter, the external magnetic 
field determines the boundary conditions of the problem. As a result, the answer to the 
question of completeness of Meissner screening depends on the results of the competition 
of contributions from a paramagnetic (spin) j and diamagnetic (charge) — 4A parts of 
the total current J. 

Similar to energy distribution in the fractional quantum Hall effect with the fiUing 
factor V — p/q and p,q G the energy gain in Eq.(19) depends on the relation |L|/|Q|. 
The Hopf invariant Q & Z numbers vacuums jHS] and is equivalent to the degree of 
degeneracity q of the ground state. The index L plays the role of the fiUing degree p of 
incomressible charged fluid state in the fractional quantum Hall effect. From this point 
of view the multiplier (1 — |L|/|(5|) in Eq.(19) is equivalent to the filling factor 1 — v for 
holes. The distinction of our system from the states in fractional quantum Hall effect 
is that the superconducting state is compressible and here (as it was mentioned) the 
effective number L of charge degrees of freedom is not an integer in the general case. 
The configurations of fields n and c = a with the integers L = Q, satisfying the relation 
of self-duality Fn — Fc, correspond to the minimum value of free energy. In this limit 

Kap is proportional to the matrix ^ | | ^ , which was used in [HEI to describe the 

topological order in the theory of fractional quantum Hall effect with the fiUing factor 
l-v at iy= 1/2. 

The boundary conditions which determine the momentum c and the topological in- 
variants L, Q depend not only on the values of the tuning parameter po of the model 
and the external magnetic field. Their physical sense and value depend also on the di- 
mensionality of the manifold for which the model is defined. In the (3 + 0)-dimensional 
case of the free energy (3) the Hopf invariant (6) is analogous to the Chern-Simons ac- 
tion (fc/47r) / dt (P X e f_,i,\a^j,dva\- This term in the action of (2 -I- l)-dimensional systems 
describes the dependence of the nonlinear modes' contribution to the free energy on the 
statistical parameter k. The coefficient k in the Chern-Simons action has the geomet- 
rical sense of the braiding number of the excitation world lines. In particular, when 
semi-fermion excitations (semions) permutate and return to the initial positions on the 
plane, the world lines braid twice and k — 2. Therewith, the statistical correlations of 
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nonlinear modes have the character of attraction and for the values k 2 their great- 
est contribution to the energy is of the order of several percents For the energy 
scale (0.1 - 1) eV this gives several tens or hundreds of degrees. Taking into account 
the relation between the dimensionality of the systems at their dynamic and statisti- 
cal descriptions, we note that the (2 + l)-dimensional case k = 2 with the open ends 
of excitation world lines is equivalent (after identification of the ends) to the compact 
statistical (3 -I- 0)-dimensional example of the Hopf linking with Q = 1. 

The (3 -I- 0)-dimensional and (2 -t- l)-dimensional situations differ by the topology of 
the regions of the field definition n |38[l3iJlf^ . When the system is periodic in one of the 
space variables and also when calculating the partition function in planar systems, one of 
three coordinates - Matsubara variable - is a periodic variable. This means that instead 
of the sphere S'^, we deal with the topology of a three-dimensional torus T'^^^ — S'^ x 
or = X X and with the corresponding mapping classes. The content of Hopf 
invariant in this case appears to be richer j^HlEniBni- For a three-dimensional torus 

the Hopf invariant is defined modulo 2q, where q is the greatest common divisor 
of the numbers {(71,(72,93} G ^. Here qi is the degree of mapping S^, where 

is the section of with the fixed i-th coordinate. Four integral numbers {qi,Q}, 
where Q is defined modulo value 2q, give us the complete homotopic classification of 
mappings with Tri[Mapq{S^ -> S'^)] = Z2, and a fixed degree q jHHl l?M liH] . 

The geometrical meaning of this modified Hopf invariant (an integer from the range 
{0, 2q — 1}) is the same. It is a linking index of the preimages of two generic points in 
S^. The cases T^xi ^^^^ rps characterized physically by different boundary 
conditions. The boundary conditions change if an angular velocity of the rotation of the 
neutral superfluid phase in He^ increases |3H1 ISHI or an external magnetic field in our 
charged system grows. Restricting the Hopf invariant change, the transition T^^^ 
promotes the appearence of the incompressible phase. 

5. Conclusion 

One can see from the given analysis that the gain of the free energy at the transition to 
the superconducting state with c 7^ arises when there is a coherent phase associated 
with the spin degrees of freedom. This phase is characterised by a pseudo-gap (5) and a 
topological order associated with linking. If the density is rather great, the momentum 
c is small and the transition to the superconducting state is not preferable. According 
to our classification, this second state is characterized by changing values of the order 
parameters p and n. The energy loss due to the term (dkp)^ may be reduced because of 
the development of one-dimensional structures. Whether these one-dimensional charge 
structures will be open, forming stripes, or closed almost one-dimensional structures in 
the form of rings, depends on the parameters of the potential 1^(^,713). In the phase 
77,3 = const with neutral spin currents the answer will depend on the value and the sign 
of the multipHer 5(^3) in the potential V{p, n^) — —bp^ + If 6 > 0, then far from Tc 
charge structures with open ends are preferable and in the case T ^ Tc we should 
prefer rings. The first experiments that verified the existence of charge structures in the 
form of rings in the underdoped phase of planar systems were described in paper j22] . 

The superconducting phase occupies on the "temperature-charge density" phase 
diagram only a part of the region, belonging to the phase with the broken antiferro- 
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magnetic order. The bounds of its existence on the phase diagram, associated with the 
characteristic values of the density p§, are determined at great po by the inequaUty (14) 
|14j and at small po these bounds depend on the critical size of the knot [23, beginning 
from which instability of the Hopf mapping arises. 

Comparing the results of this paper based on consideration of the local fields with the 
conclusions following from the BCS-like model with two species of fermions, we pay 
attention to the following qualitative coincidence. One can conclude from Eq.(3) that the 
parameter po determines the value of the coupling constant. Therefore, the appearance 
of the solutions (different from the standard BCS model) for the superconducting gap in 
the paper with the finite value of the coupling constant is analogous to the existence 
of the threshold for small values of po in this paper. 

As distinct from the model [E] , the states considered in the present paper are signifi- 
cantly inhomogeneous. The analysis of the state Fik 7^ 0, p 7^ const is an open problem at 
present. Here we mention only that the superconducting current with the amplitude cq, 
flowing round the rings (8), gives the additional term CqR"^ to the multiplier in Eq.(lO). 
This explains in particular why the superconducting region on the "temperature-the 
level of doping" phase diagram is shifted to the Hue 5T{x) = of the transition to the 
state with the spin pseudo-gap. Indeed, in this case VefjipjU^) — —bejfP^ + {d/2)p'^ 
with beff = b — [no + Cq) — {const / B?)5T . Therefore the finite value of the momentum 
c of superconducting pairs decreases 5T. Besides, the contribution to the free energy in 
the inhomogeneous state due to {dkp)^ ^ decreases the gain in Eq.(19). 

The superconductivity in cluster systems may be evidently studied, basing on the ap- 
proaches beyond the mean field theory. For example, we may use the exact Richardson 
solution and Bethe ansatz equation jll], as well as the methods of the conformal field 
theory 02]. The exact solution of the ground state problem under the condition of the 
finite value of the total momentum of pairs at such an approach is one of the important 
problems. Since the conformal nature of the dimensionality 3/4 ^QI in (5) and (14) in- 
fluences the character of the scale which enters into energy dependent response functions 
(which is proportional to T^/"^ ESDi it should also be studied carefully. 

We are grateful to A.G. Abanov, S.A. Brazovsky, S. Davis, L.D. Faddeev, S.M. 
Girvin, V.E. Kravtsov, E.A. Kuznetsov, B. Lake, A.I. Larkin, A.G. Litvak, C. Renner, H. 
Takagi, V.A. Verbus, G.E. Volovik, Yu Lu, Y.-S. Wu for advices and useful discussions. 
This paper was supported in part by the Grant of the Russian Foundation for Basic 
Research No. 01-02-17225. 



Appendix 

The proof of the Eq. (14) uses the following chain of inequalities: 

|L|<||C||6 . ||H||6/5^61/S||[VXC]||2 • ||H||6/5^ (20) 

6V6||[vxc]|i2 • ||H||?/^ . wml''^ 

(327r2)-4/3FV2i;^2/3^y6 ^ [2,2^^)-^'^ fI'^ F^/^ . 

Here ||H||p = (/ (Px\H.\pY'^ At the first and third steps we used Holder inequality: 
||f • g|| < ||f||p • ||g||g with 1/p -I- \/q = 1. Under the condition V • c = we employs 
at the second step the Ladyzhenskaya inequality 04] [40]: ||c||6 ^ 6^/^||[V x c]||2. The 
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fourth step in the set of inequahties arises after the comparison of the terms 1 1 [V x c] 1 1 and 
||H|| with the terms Fn and Fc in Eq.(13). The last Hne shows also separate contributions 
from n- and c-parts of the free energy (13) to the finite result (14). Using a chain of 
Holder and Ladyzhenskaya inequalities, analogously one may find that: 

J^y'^c^' ^ (167r2)4/3|L'|. 

The coefficient in this inequality differs from (14) due to the coefficient 1/4 (because of 
the charge 2e of superconducting pairs) of the first term of the free energy Fc in (13). 
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